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ABSTRACT 

Context. The census of Galactic globular clusters (GCs) is still incomplete, and about ten new objects are supposed to await discovery, 
hidden behind the crowded and reddened regions of the Galactic bulge and disk. 

Aims. We investigated the nature of three new GC candidates, discovered in the frames collected by the Vista Variables in the Via 
Lactea (VVV) near-infrared survey. They will be called VVV CL002, VVV CL003, and VVV CL004. 

Methods. We studied the results of point- spread-function near-infrared photometry from VVV data for the three objects and their 
surrounding fields, the proper motion information available in the literature and, when possible, we derived the cluster parameters by 
means of calibrated indices measured on the color-magnitude diagrams. 

Results. The evidence shows that VVV CL002 is a newly discovered, small, moderately metal-rich ([Fe/H] « -0.4) Galactic GC. It 
is located at a Galactocentric distance of 0.7+0.9 kpc, and it could be one of the nearest GC to the Galactic center. Its characteristics 
are more similar to those of low-mass, Palomar-like GCs than to more classical, old, and massive bulge GCs. VVV CL003 is the 
first star cluster discovered in the Galactic disk on the opposite side of the center with respect to the Sun, at a Galactocentric distance 
of ~5 kpc. Its high metallicity ([Fe/H] « -0.1) and location point to an open cluster, but a GC cannot be excluded. VVV CL004, 
on the contrary, is most probably only a random clump of field stars, as indicated by both its low statistical significance and by the 
impossibility to distinguish its stars from the surrounding field population. 

Conclusions. We claim the detection of i) a new Galactic GC, deriving an estimate of its basic parameters; ii) a stellar aggregate, 
probably an open cluster, in the disk directly beyond the Galactic center; and Hi) an overdensity of stars, most probably an asterism. 
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1. Introduction 

The census of Galactic globular clusters (GCs) is still incom- 
plete, because distant objects can easily hide behind the very 
crowded and highly reddened stellar fields in the direction of the 
Galactic bulge and disk. Indeed, [Ivanov et aL] ( |2005| ) estimated 
that there should be 10+3 undiscovered GCs hiding in the inner 
Milky Way. The ad vent of the new generat ion of exte nsive sur- 
veys such a s SDSS ( [Abazajian et al.|2009>, 2MASS ( |Skrutskie| 
|etaL||2006) , and GLIMPSE ( [Benjamin et al^|2003l ) permitted 
the detection of sev eral new Galac tic GCs. The December 2010 
compilation of the Harris (1996 ) catalog included seven new 
GCs not present in the February 2003 version, but seven more 
objects have been proposed in th e last years: SDSSJ1257+3419 
( [Sakamoto & Hasegawa| [2006]) FSR 584 pica et al.||2007]>, 
FSR17 67 (Bonatto et al.||2007), FSR 190 { proebrich et al. 



2008a), Pfl eiderer 2 ( |Ortolani et ai.|2009), VVV C L001 ([Minniti 



et al.|201l) , and Mercer 5 ( Longmore et al.|20TT) . The family of 



Galactic GCs is thus steadily increasing, although the true nature 
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of the GC candidates is not always easy to unveil, and some ob - 
jects are still debated (see, for example, [Froebrich et al.|2 008b ). 

Globular clusters are old objects, whose age ranges from 
about half to nearly one Hubble time. Their spatial distribution 
and physical/chemical properties thus bear information about the 
Galactic formation and evolution processes. However, they con- 
tinuously loose mass during their entire lifetime because of both 
external processes such as Galacti c tidal stress, and internal ones 
such as stellar evaporation (Meylan & Heggie|1997[ |Gnedin & 



|Ostriker|[l997| ). The identification and the study of low-mass 
GCs in the Galactic bulge can therefore provide new informa- 
tion about both its formation and chemical enrichment history, 
and the tidal forces governing its dense environment. 

In this paper we analyze three new cluster candidates dis- 
covered in the VISTA Variables in the Via Lactea (VVV) Public 
Survey, wh ich is gathering near-IR data of the Galactic bulge and 
inner disk ( [Minniti et aLl|2010[ |Saito et al.||2010[ |Catelan et aL] 
|2011| ). This survey is scanning the central regions of the Milky 
Way with unprecedented depth and extension, thus generating 
the ideal database to search for missing GC s in the inner Galaxy . 
Following the nomenclature introduced by Minniti e t al.| ( [20 1 1 1 ) 
and Boriss ova et al.| ( |20TT) , the objects under study here will be 
called VVV CL002, VVV CL003, and VVV CL004. 
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2. Photometric data 

2.1. The cluster candidates 

VVV CL002, VVV CL003, and VVV CL004 were identified 
as previously unknown stellar overdensities by visual inspection 
of the three-band (/, H, K s ), color composite frames produced 
by the VVV survey. Table [T] reports for each object the corre- 
sponding VVV fi elds and the J2000 coordinates of their center, 
as determined in § |3.1| Figure [T] shows the three objects in the K s - 
band stacked tiles downloaded from the Vista Science Archive 
websitcQ(VSA). VVV CL002 and VVV CL003 look like loose 
and low-mass, Palomar-like GCs, although VVV CL003 is less 
clear because it is embedded in a very crowded stellar field. 
VVV CL004, on the contrary, does not show a clear concen- 
tration of faint stars, but only a local overdensity of a few bright 
objects. 
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Fig. 1. 5' x 5' logarithmic gray-scale K s -bsnid images from 
the VVV archive, centered on the three candidate clusters 
(from upper to lower panel: VVV CL002, VVV CL003, and 
VVV CL004). North is up and East is to the right. The white 
circles, centered on the derived cluster center, indicate the half- 
light radius estimated for each object. 



2.2. Data reduction and photometry 

After the discovery, we retrieved from VSA the stacked im- 
ages of the individual 2048x2048 pixel exposures containing 
the three candidate clusters. The quantity of frames used in each 
band is indicated in Table [T] The effective exposure time of the 
H- and ^ s -band images was 8s, while it was 24s in J. The frames 
in the Y and Z bands, included in the observation plan of the 
VVV survey, were not available at the time of the investiga- 
tion. The data were acquired on 2010 April 7 to 14, with the 
VIRCAM camera m ounted on the VISTA 4m telescope at the 
Paranal Observatory ( [Emerson & S utherlan d|2010| ), and reduced 
at the Cambridge Astro nomical Survey U nit (CASUjjwith the 
VIRCAM pipeline vl.O ( [Irwin et al.|2004| ). During the observa- 
tions the weather conditions fell within the survey's constraints 
for seeing, airmass, and Moon distance (Minniti et al.|2010"| ), and 
the quality of the data was satisfactory. 

Stellar photometry was performed with the DAOPHOTII 
and ALLFRAME codes ( |Stetson|[T994l ). After the construction 
of the point-spread function (PSF) for each frame and an initial 
profile-fitting photometry, all images of the same object were 
matched to create a combined frame, where the master list of 
star-like sources was created. A PSF fitting of this master list was 
then performed on each single exposure. The instrumental mag- 
nitudes were transformed into the 2MASS photometric system 
by means of the 2MASS stars found in each field. For each frame 
the 2MASS sources were first identified in our catalog, then the 
list was cleaned of objects that in our photometry had a detected 
source nearer than 2'/2 and fainter by less than 3.88 magnitudes. 
We thus excluded the stars whose 2MASS aperture photometry 
could have been contaminated by nearby undetected sources by 
more than 0.03 magnitudes. We also limited the selection of cal- 
ibration stars in the range /2mass=1 1.5-13, //2mass= 10-11, and 
^2mass,s=9.5-11, to avoid stars heavily saturated in VVV data 
and the fainter end of the 2MASS survey, where systematic er- 
rors can be significant (see below). Between 70 and 490 stars 
met all criteria in each frame. The calibration equations were in 
the form ra 2 MASS - ra=A+Bx(/- J K s ) 2 MASS, where ra 2 MASS and m 
are the 2MASS and instrumental magnitudes, respectively, and 
A and B are constants. The color term is in general negligible, 
being on the order of few hundredths of magnitude. In Table [2] 
we give the average coefficients A and B for each pair of offsets 
collected within a few minutes of each other, and the rms of the 
points around the fitted relation. The final calibrated photometry 



1 http://horus.roe.ac.uk/vsa/ 

2 http://casu.ast.cam.ac.uk/ 
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Table 1. Coordinates of the center of the three GC candidates, and the VVV data used in the present investigation. The classification 
(GC= globular cluster, OC= open cluster, NC= not a cluster) resulting from our investigation is given in the last column. 



ID 


RA 


dec 


/ 


b 


VVV field 


frames 


classification 




hh:mm:ss 


deg:mm:ss 


deg 


deg 








VVV CL002 


17:41:06.30 


-28:50:42.3 


359.5586 


0.8888 


b347 


2x7, 2x//, 4xK s 


GC 


VVV CL003 


17:38:54.56 


-29:54:25.3 


358.4047 


0.7298 


b346 


2x7, 2x7/, 4xK s 


OC? 


VVV CL004 


17:54:31.91 


-22:13:37.8 


6.7900 


1.7179 


b352 


4x7, 4x7/, 4x7: s 


NC 




Fig. 2. Magnitude difference between 2MASS and our photom- Fig> 3< Comparison between our 7-band magnitude (7 our ) and the 

etry for all the stars in common as a function of 2MASS magni- V SA aperture photometry catalog (7 VS a), for stars in the stellar 

tude in the stellar field of VVV CL002. From top to bottom: 7, fie i d of V VV CL002. The gray line indicates the zero-difference 

H, and K s magnitudes. between the two photometries (/vsa=/out). 



was checked studying the magnitude difference with respect to 
the 2MASS catalog. These residuals showed no trend with color, 
but a systematic difference was found for the fainter stars in all 
bands and fields. As an example, the results of this comparison 
for VVV CL002 and its stellar field is shown in Figure [2] but 
very similar plots were obtained for the other two cluster candi- 
dates as well. The numerous and dispersed negative values indi- 
cate the 2MASS stars contaminated by undetected sources. The 
agreement is satisfactory for the brighter stars, but in the faintest 
~1.5 magnitudes the bulk of the distribution curves upward to 
positive values, indicating that the stars are fainter in the 2MASS 
catalog than in our photometry. This feature was found identical 
even when performi ng PSF photome try of VVV data with the 
DoPhot code ( [Schechter et al.|[T993] ), and when comparing the 
2MASS magnitudes with the VVV aperture photometry down- 
loaded from VSA, whose catalogs agree well with our results. As 
an example, in Figure [3] our 7 magnitudes in the VVV CL002 
stellar field are compared to the corresponding quantity in the 
VSA catalog. The same data that show a systematic difference 
when compared to 2MASS (upper panel of Figure [2]) do not re- 
veal any trend in this plot, apart from an expected mismatch 
for stars in the saturation regime (7 < 12.5), where the VSA 
aperture photometry underestimates the stellar brightness. We 
conclude that this systematic offset reflects an intrinsic differ- 



ence between 2MASS and VVV data. We believe it is caused by 
the relative shallowness of the former and its poorer pixel scale, 
which leads to an overestimate of the sky background, especially 
in crowded fields, and consequently to an overly-faint magnitude 
estimate for faint stars. A more detailed analysis is beyond the 
scope of this paper, and it will be presented in a future study 
(Mauro et al. 2011, in preparation). 

VVV CL003 was found in frames acquired with the chips 
#12 and #16 of the VIRCAM camera. The latter chip is never 
correctly flat-fielded because of short-timescale variations of its 
quantum efficienc)^] We found that the PSF photometry result- 
ing from this chip was unreliable, because of large zero-point 
variations on the CCD. We finally decided to exclude these 
frames, and only the exposures collected with chip #12, listed 
in Table [T] were used. Unfortunately, the object fell on the edge 
of this detector, where the photometric results are of lower qual- 
ity, and the cluster center was only ~ 19" from the border. Only 
half of the surrounding field was therefore available for the anal- 
ysis, and we estimated by eye that we have studied only about 
three-fourths of the candidate cluster. 



3 http://casu.ast.cam.ac.uk/surveys-projects/vista/technical/known- 
issues 
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3. Results 

3.1. Cluster center 

With the final photometry of the three cluster candidates, we re- 
vised the location of their centers, first estimated by eye. This 
calculation was affected by a smooth background of numerous 
faint field sources that shifted the barycenter of the stellar dis- 
tribution toward the center of the region selected for the calcu- 
lation. Cuts in the color-magnitude diagrams did not help, be- 
cause either the distribution of cluster and field stars was sim- 
ilar (VVV CL002 and VVV CL004, see Figure §>, or the pho- 
tometric information was not reliable (VVV CL003, see g2.2| ). 
However, if r is the position of the center of the stellar distribu- 
tion, and N is the number of stars, we can write 



Table 2. Average calibration coefficients for each pair of con- 
tiguous offsets. 



r*N* = ^n = r c N c + r f N f , 



(1) 



where r, is the position of the z-th star, and the subscripts *, f, and 
c refer to all, field, and cluster stars, respectively. After simple 
algebra, considering that TV* = N c + Nf, and that the center of 
the distribution of field stars (r f ) coincides with the center of the 
area included in the calculation (r ), the previous relation yields 

Z/(r/ - r ) 



r c = r + 



N*-Nf 



(2) 



This formula was used to find the position of the centers of the 
three objects. The total number of field stars (Nf) was estimated 
from the average stellar density in regions farther than 3' from 
the initial by-eye center position. To limit the influence of ran- 
dom fluctuations of the stellar distribution, the calculation was 
repeated in 360 areas of 3' radius, whose centers were equally 
spaced on a circle of 0f5 radius around the initial guess. The re- 
sults differed by less than \" with an rms of only Off 3, and they 
were averaged to obtain the final cluster center, given in Table [T] 
To cover a larger field for VVV CL003, the stellar counts 
were performed including the data collected with chip #16. 
As already noted, these data are not suitable for photometry - 
though they are still good enough for star count purposes. 

3.2. Significance of the stellar overdensity 

Following Kop osov et al.| ( |2007] ), we can estimate the statistical 
significance of the stellar overdensities from the number of stars 
detected in excess to the local background, whose random fluc- 
tuations are assumed to be Poissonian. We detected 6965 stars in 
the inner circle of radius V of VVV CL002, while the field den- 
sity level can account only for 5744 objects: an overdensity of 
stars at the 10.9<x level is thus present in this area. Analogously, 
5193 stars are detected for VVV CL003, against a field expec- 
tation of 4313, and the overdensity is significant at the 9cr level. 
The quantity of stars detected in excess of the field expectation 
within V from the center of VVV CL004 is lower (-400), and 
the statistical significance is ~4.5<x. The 520 square degrees cov- 
ered by the VVV survey can be divided into -600, 000 circles of 
radius V, therefore the probability of detecting one or more such 
overdensity in the whole survey that is caused by random field 
fluctuations is -98%, and about four similar objects should be 
expected. Hence, the statistical significance of VVV CL004 is 
low, casting serious doubts on its nature as a true stellar cluster. 

3.3. Radial density profile 

After the determination of the cluster center, we calculated the 
stellar density (number of stars per arcmin 2 ) in concentric rings 



date 


UT 


filter A 


B 


rms 






mag 




mag 


VVV CL002 - chip 12 


2010-04-07 


7:17:28 


H 0.17+0.03 


-0.013+0.011 


0.057 


2010-04-07 


7:20:46 


K s -0.47+0.02 


-0.060+0.009 


0.061 


2010-04-07 


7:24:40 


J 0.64+0.02 


0.009+0.008 


0.053 


2010-04-14 


9:32:23 


K s -0.56+0.02 


-0.028+0.008 


0.052 


VVV CL003 - chip 12 


2010-04-07 


6:52:54 


H 0.23+0.02 


-0.014+0.010 


0.053 


2010-04-07 


6:56:11 


K s -0.561+0.016 


-0.017+0.007 


0.058 


2010-04-07 


7:01:05 


J 0.55+0.02 


0.039+0.008 


0.041 


2010-04-14 


9:06:35 


^ s -0.559+0.014 


-0.019+0.006 


0.053 


VVV CL004 - chip 5 


2010-04-11 


5:41:06 


H 0.278+0.015 


-0.048+0.010 


0.033 


2010-04-11 


5:44:13 


K s -0.554+0.009 


-0.025+0.006 


0.036 


2010-04-11 


5:48:40 


J 0.592+0.011 


0.039+0.008 


0.032 


VVV CL004 - chip 9 


2010-04-11 


5:40:12 


H 0.230+0.010 


-0.014+0.007 


0.028 


2010-04-11 


5:43:20 


K s -0.563+0.010 


-0.014+0.007 


0.032 


2010-04-11 


5:47:07 


J 0.622+0.011 


0.020+0.007 


0.030 



between the circles of radius Of 05 x N and Of 05 x TV + Of 1 (with 
N=0, 1,2,3...). The resulting radial profiles of the stellar density 
are shown in Figure [4] where the error bars are given by the 
Poissonian noise. The field density level, indicated by a dotted 
line, was assumed as the average value between If 8 and 2' from 
the center, i.e. beyond the clusters' tidal radii determined below. 
All three objects show a clustering of stars in the central region. 
The innermost bin shows a sudden increase of the density in all 
cases, but this value is not very reliable because it corresponds 
to a circle of only Of 1 radius. On the other hand, the possibility 
of photometric incompleteness in the most crowded region must 
be taken into account in the following analysis. If the catalogs 
were incomplete in the inner part of the cluster candidates, the 
statistical significance of the detected stellar overdensity would 
be underestimated, and the decontaminated color-magnitude di- 
agrams would be relatively deficient of stars, mainly at fainter 
magnitudes. However, we detect between 80 and 100 stars in the 
inner Of 1 of the three clusters, that is 225 times the area covered 
by a circle of diameter equal to the FWHM of the PSF (078). 
This indicates that the three objects are not excessively crowded, 
and the incompleteness should not be very severe. 

The radial profiles of VVV CL002 and VVV CL003 are very 
similar, smoothly falling off at increasing distance into an ex- 
tended tail. On the contrary, the density of VVV CL004 drops 
more abruptly, but the cluster limit is not immediately evident 
because of strong density fluctuations. All these conclusions are 
confirmed and better visualized by the radial profiles of the cu- 
mulative fraction of stars, presented in Figure [5] These were cal- 
culated counting the stars detected within a given radius, and 
subtracting the number of expected field stars, estimated from 
the average density between If 8 and 2' from the center. This 
range was also used to derive the total number of cluster stars 
that were used to normalize the curve. The resulting profiles 
show that VVV CL002 extends out to ~ lf6 from the center, 
and its half-light radius, i.e. the radius containing half of the de- 
tected stars, is rh ~ Of 75. The apparent size of VVV CL003 
is probably smaller because it fades into the field density be- 
tween lf4 and lf6 from the center, and its half-light radius is 
^ ~ Of 6. The angular size of VVV CL004 is more uncertain, 
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because the field density is patchy and both the density profile 
and the cumulative fraction of stars fluctuate strongly. The ra- 
dius where the density merges into the field continuum is un- 
clear, probably comprised between Y and If 7, while r h ~ Of 35. 
The radial d ensity profile of the three objects was fitted with a 
King| < |1966| > profile, and the results are overplotted on the ob- 



servational data in Figure [4] For VVV CL003 this procedure re- 
turned a very good fit of the data, with a tidal radius r t = If 8 
and a central concentration c = log (r t /r c ) = 0.56, where r c is 
the core radius. The density profile of VVV CL004 is also sat- 
isfactorily reproduced by the fitted function, and the cluster pa- 
rameters derived in this way are r t = If 7 and c = 0.75. On the 
other hand, the routine failed to return physically reasonable pa- 
rameters for VVV CL002, because the density between Of 5 and 
If 5 is too high, causing r t to diverge. A shallowe r slope of the 
density profile indicates the presence of tidal tails ( Kupperetal 



|2010b| ), and the cluster may be losing stars along its orbit. We 
therefore estimated r t = If 8 from the point where N(r<d)/N tot =l 
in Figure [5] then a King profile of fixed tidal radius w as fitted 
to find c = 0.65. Similar problems were found with a|Wilson 
( |1975| ) profile, calc ulated from the ene rgy distribution function 
through the code of Sollima et al. (2009 ): fitting both the central 
decay (r < 0f5) and the enhanced tail was impossible. In this 
case the fitting routine converged, but to a very high tidal radius 
(c=0.68, r t = 4f8), and the r esulting fit was u nsatisfactory. We 
found that a EFF power-law (Elso n et al.|1987| ) fits the observed 
profile well, with central density //o=5630 stars, normalization 
radius a = Off 16, and power -0.086 (y=0.17). The fitting rou- 
tine likely underestimates the uncertainties associated with its 
output, because it does not take into account any source of sys- 
tematic error, mainly the definition of the background level. We 
consider that the width of the bins of the density profile (Of 1) is 
a better estimate of the errors on the core, half-light, and tidal 
radii. Indeed, we found that these quantities can vary by up to 
Of 1 from the given values when different but still reasonable as- 
sumptions on the background level were adopted. 



3.4. Color-magnitude diagrams 

The color-magnitude diagrams (CMDs) of the three objects are 
shown in the left panels of Figure |6| where the cluster stars are 
overplotted on the photometry of a nearby comparison stellar 
field. The cluster area of VVV CL002 and VVV CL004 was 
limited to a circle of radius equal to the half-light radius to min- 
imize the background contamination. The comparison field was 
an annulus of the same center and area as the cluster, and inner 
radius of If 8, i.e. just outside the cluster tidal radius. A sim- 
ilar definition of cluster and field areas for VVV CL003 was 
complicated because its photometric catalog lacked the required 
symmetry about the center. We adopted as cluster region a cir- 
cle of radius 0f5, but part of this area fell outside the CCD. The 
definition of the comparison field was also complicated because 
VVV CL003 is found in a region affected by a very strong red- 
dening gradient. Indeed, the interstellar extinction varies by up 
to AE(B - V) «2 over the 10' of the chip, as confirmed by in- 

(|1998|) maps. We noticed, how- 



spection of the |Schlegel et al, 



ever, that the CMD of the cluster area presents two well-defined 
red giant branches (RGBs): a narrow and well-populated one at 
(/ - K s ) «2, and a broader RGB about 0.5 magnitudes redder. 
The bluer sequence was found in all studied areas, always very 
similar in terms of quantity of stars, width, and magnitude of the 
red clump, although shifted to different colors. The other RGB 
component, on the other hand, was not observed outside ~ Y 
from the cluster center. We therefore concluded that the bluer 



4000 



3000 



2000 - 



£2000 - 



3000 



2000 




d (arcmin) 

Fig. 4. Radial profile of the stellar surface density for the three 
cluster candidates. The dotted lines show the estimated field 
level, while the curve indicates the best-fit King profile. 



VVV CL002 



0.5 - 



0.5 




3 



d (arcmin) 

Fig. 5. Cumulative fraction of cluster stars as a function of dis- 
tance from the center. The field and the half-light levels are indi- 
cated with full and dotted lines, respectively. 



RGB represents the field population, and we searched for a re- 
gion in the vicinity of the cluster where the color of this feature 
coincided with that observed in the cluster area. The adopted re- 
gion had the same shape and dimension as the defined cluster 
area, but offset Y toward the northwest, i.e. perpendicularly to 
the Galactic plane toward increasing Galactic latitude. Part of 
the object could still fall in this area, but its density should be 
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very low at this distance. In Figure[7]we show the color-color di- 
agrams of the three objects and their comparison fields. This plot 
was restricted to stars brighter than K s =15 to highlight the color 
differences at the RGB level, because the increased photometric 
errors blur them for fainter stars. 

As already noted, the CMD of VVV CL003 shows two dis- 
tinct RGB sequences, and the redder one is not observed in the 
surrounding stellar field. In the comparison field, some residual 
cluster stars are also visible, not observed at longer distances 
from its center. On the other hand, the CMDs of VVV CL002 
and VVV CL004 show no immediately clear difference with 
their comparison field, except for the enhanced number of stars. 
Further inspection reveals that the cluster RGB of VVV CL002 
is on average bluer than the field population. This is visible even 
in the color-color diagrams, where cluster stars of VVV CL002 
are offset to slightly bluer colors with respect to the field popula- 
tion. We found no color gradient that can explain this behavior in 
the whole ~ 10' x 10' area covered by our photometry, although 
the reddening pattern is probably variable on a small scale, as in- 
dicated by some highly reddened stars found in both cluster and 
comparison fields. VVV CL004 is less clear than the other two 
objects. Its color-color diagram contains more red stars than the 
surrounding field, but this is a simple consequence of the higher 
number of stars in the upper (redder) RGB of the cluster area. 
Indeed, although some bright (K s = 11-13) red stars without field 
counterparts are observed, the CMD of VVV CL004 shows no 
clear peculiarity with respect to the field, except for being much 
richer in brighter stars, as well as in faint red stars. 

3.5. Statistical field decontamination 

The CMDs of Figure [6] are strongly contaminated by field stars, 
and so we applied a statistical procedure to isolate the clus- 
ter CMDs. The code written f or this purpose was based on the 
method of Gall art et aT] ( |2QQ3| ). In brief, a metric was defined in 



the K S -(J - ^ s ) plane as 



^(K s , a -K s , h ) 2 + (kx((J-K s ) a -(J 



(3) 



where k is an arbitrary coefficient weighting a difference in color 
with respect to a difference in K s . The code scanned the list of 
stars in the comparison field and found the nearest star (smaller 
d) in the CMD of the cluster area for each object. This source 
was rejected as a contaminating field star if its distance d was 
smaller than an arbitrary threshold d max , otherwise the star in the 
comparison area was flagged as an "object without a counter- 
part". In the first case, the list of stars in the cluster area was 
updated after the rejection before proceeding to the next entry in 
the list. 

The choice of th e parameters d max an d k is not completely 
free of arbitrariness. |Gallart et al.| ( |2003| ) used d max =03 mag- 
nitudes, and we found this is a good choice in our case also, 
because it prevented the association between two stars with a 
very different position in the CMD. The same authors adopted 
k-1 in their optical photometry, which gives a dominant role 
to an offset in color. Their choice was dictated by the need to 
remove foreground contaminants distributed in a wide range of 
distances, with only limited variations of reddening. In our case, 
though, the patchy and locally variable structure of the reddening 
makes a color offset between cluster area and comparison field 
more likely than a difference in magnitude. Moreover, (J - ^ s ) 
is more sensitive than the K s band to interstellar absorption, and 
even a different distance can affect the stellar color in our case, 
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because we are observing through the dense clouds of the in- 
ner Galaxy. In conclusion, a much lower value of k should be 
appropriate in our study. Indeed, we explored the effects of this 
choice through repeated trials, searching for the value that re- 
turned a clearer CMD and a smaller fraction of comparison field 
stars that were without counterpart. We found that both the resid- 
ual contamination in the product CMDs and the quantity of un- 
matched field stars decreased with k, and k <2 was clearly pre- 
ferred. Only negligible changes were observed when the scaling 
factor was varied in the range k=l-2, and we finally adopted 
£=1.6. The unmatched field stars thus decreased to ~2% for 
VVV CL004, which is below the Poissonian noise, while it was 
-10% for VVV CL002 and VVV CL003, i.e. about 3^ times 
the Poissonian noise. However, the decontaminated CMDs of 
VVV CL002 and VVV CL003 contain 453 and 211 stars, re- 
spectively, which matches the expectation of 454 and 250 ob- 



jects well, that was a result of the stellar counts of Section 3.3 



Hence, the decontamination procedure worked satisfactorily. 

In the decontamination procedure, different choices of the 
cluster and field areas necessarily led to different sets of removed 
stars. However, all features observed in the resulting CMDs, 
such as sequences and clumps, were observed in all cases when 
the cluster center was unchanged. The innermost regions of the 
three objects were excluded from the calculation, because our 
photometric catalogs could be incomplete in the most crowded 
regions. In general, the cleaning procedure left some residual 
field contamination, that decreased when choosing i) smaller ar- 
eas, which consequently reduced the number of contaminating 
stars to be removed, and if) a comparison field closer to the clus- 
ter center, which minimized the effects of a different local red- 
dening. In both cases, the lower level of residual contamination 
came at the cost of a poorer derived CMD, because cluster mem- 
bers were also excluded. Nevertheless, these choices are justi- 
fied here, because our aim is to unveil the behavior of the clus- 



ter CMD, while the true number of stars in the different CMD 
branches (i.e. the detailed luminosity function) is of lesser im- 
portance. Finally, the cluster area adopted in the decontamina- 
tion procedure of VVV CL002 and VVV CL004 was a circle of 
radius Of 5 and Of 3, respectively, slightly smaller than their half- 
light radii. The inner Of 1 from the center was also excluded. The 
comparison fields were annuli of equal area as the cluster region, 
with inner radius of 0f9 in the case of VVV CL002 and Of 7 for 
VVV CL004. The cluster density is very low at these distances, 
and the loss of cluster stars in the resulting CMDs should be lim- 
ited. Yet the a reas of VVV CL003 and its comparison field were 
defined as in g3.4| (two circles of Of 5 radius, offset by 1') because 
the strong underlying reddening gradient and the position of the 
cluster at the edge of the CCD gave no freedom for exploring 
different choices. Even for VVV CL003 the inner region with 
r < Of 1 from the center was excluded. 

The decontaminated CMDs of the three cluster candidates 
are shown in the right panels of Figure [6] In all cases a clear 
evolutionary sequence of stars is visible, with a clumpy behavior 
caused by some underpopulated regions, and many highly scat- 
tered stars mainly at the fainter end of the diagrams. These fea- 
tures are probably an artifact of the decontamination procedure: 
when the cluster and field CMDs are very similar and the sam- 
ple is statistically small, stars in crowded regions of the CMD 
are more easily matched with a field counterpart than the ob- 
jects deviating because of high random error or higher differ- 
ential reddening. As a consequence, the algorithm could over- 
subtract stars along the stellar sequence and under- subtract them 
in the wings of the color distribution, while still subtracting the 
expected amount of field contaminants, as indicated by the quan- 
tity of cluster and unmatched field stars given before. 
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3.6. Reddening-free indices 

In the presence of non-uniform, spatially variable reddening, it 
is instructive to investigate the CMDs where the color and mag- 
nitude are substituted by reddening-free indices. In Figure [8] we 
show the decontaminated CMDs of the three objects overplot- 
ted on their comparison field, with the indices C3 = (J - H) - 
1.47 x(H- K s ) and m 4 = K s - 1.22 X (J - H) in pla ce of the 
color and m agnitude. Th ese quantities were defined by Catelan 
etal.| ( |20lT| assuming the Cardell fet al.| ( [1989] ) relation b etween 
extinction in different bands, and R=A V /E(Z? - V)=3.09 (Rieke 



& Lebofsky||1985| ). We also experimented with analogous in- 



dices derived assumi ng the near-IR extint ion law for the cen- 
tral Galactic regions ( [Nishiyama et al.|2009| ), but they led to the 
same general conclusions. Unfortunately, the error on C3 and 
blurs the resulting diagram, because it is about three times larger 
than that in (J - K s ) and K s respectively, owing to propagation 
of errors. The average uncertainty on these indices is about 0.1 
mag down to ra4=10, 0.15 mag at ra4=12, and it increases ex- 
ponentially at fainter magnitudes. Thus most of the information 
is erased. Moreover, while zeroing the effects of reddening, the 
indices also strongly limit the effect of other parameters, such 
as age or metallicity. For example, we verified through synthetic 
isochrone experiments that the RGBs of two stellar populations 
of the same age but different metallicity, with [Fe/H]=0.0 and 
-1.5, would form two nearly parallel sequences separated by 
less than AC3 »0.2. This, combined with the large errors, makes 
two different populations difficult to distinguish. Nevertheless, 
Figure [8] shows that the cluster sequence is redder than the field 
in VVV CL003. For VVV CL002 the difference is less evident, 
but both a linear and a quadratic fit of the upper RGB {m^ < 12.5) 
are between 0.01 and 0.04 magnitudes bluer in the cluster than in 
the field. Then, the color differences between these objects and 
their surrounding fields are not caused by differential reddening 
alone. 

Figure [8] does not reveal much about VVV CL004, because 
field and cluster sequences already overlapped in the (J-K s )-K s 
plane. It can be noted, however, that some red stars in the upper 
RGB (K s = 11 - 13, see Figure [6} agree with the field locus in 
the C3-m4 diagram, indicating that some differential reddening is 
present even in the direction of this object. 



3.7. Proper motions 

The information about proper motion of stars in the field of an 
observed stellar overdensity is a very powerful tool to clarify its 
nature, especially if combined wi th the photometric data (e.g., 
Carraro 2005 ; Pavani & Bica 2007 ). Indeed, a study of new clus- 
ter candidates without an analysis of the proper motions would 
seem incomplete. Unfortunately, it provides no real information 
for our three objects, because available data are sc arce, and at 
the estimated distances (~10 kpc, see g3.9| - g3.10] ) a too high 
precision is required to detected an appreciable diff erence with 
the fie ld population. For example, the UC AC3 (Zacharias et al.| 
20 10t and NOMAD ( [Zacharias etal.|2004| ) catalogs contain only 
2 and 1 1 stars, respectively, in the inner V of VVV CL002. 

The proper motions of stars in the direction of VVV CL003 
and VVV CL004 and in their surrounding fields drawn from the 
PPMXL catalog (Roeser ~et al.||2010| ) are compared in the mid- 
dle and lower panel of Figure |9| The cluster stars were taken 
within 0f8 (VVV CL003) and 0f5 (VVV CL004) from the cen- 
ter, and the comparison field was defined as an annulus between 
2' and 5'. Clearly, this comparison is inconclusive, and the av- 
erage value in the cluster and field area do not differ by more 
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Fig. 9. Distribution of proper motions in the cluster area (full 
dots with error bars) and in the comparison field (small gray 
dots). Upper panel. VVV CL002; Upper panel. VVV CL003; 
Lower panel. VVV CL004. 
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than 1 mas yr _1 in both components. According to our stellar 
counts, only 11% and 13% of the stars in the selected region 
should be members of VVV CL003 and VVV CL004, respec- 
tively. Attempting to isolate cluster members of VVV CL003 
through the 2MASS photometry given in the PPMXL catalog, 
only three stars compatible with the cluster sequence can be 
identified. 

The proper motions of stars within Of 8 from the center of 
VVV CL002 are offset with respect to the comparison field, as 
shown in the upper panel of Figure [9] Only 4.4 cluster mem- 
bers (~20%) are expected among the 22 entries of the PPMXL 
catalog, and its very hard to assess the significance of the ob- 
served behavior. As an exercise, we will check if this offset 
is compatible with the detection of the peculiar movement of 
VVV CL002. Assuming the escape ve locity from the Galactic 



bulge (~ 600 km s~\ Lamb et al. 1996| ) as an upper limit for 



the tangential velocity, the proper motion of stars at 7.5 kpc 
from the Sun should not exceed 17 mas yr _1 . We will there- 
fore restrict the calculation to the stars with <20 mas yr _1 
because, considering the average error (~8 mas yr _1 ), the clus- 
ter stars should not be found outside this range. The average 
values are QI a ,/Is)=(-2.8,-0.5) and (0.8,-3.0) for the cluster 
and field stars, respectively. In the cluster area we can write 
N tot j2 = Nfflf + N c fl c for both components of proper motion, 
where the subscripts / and c denote the field and cluster mem- 
bers, respectively, and N is the number of stars. We have N tot =12, 
and assuming iV c =5, Nf=7, we obtain (jJ a , fl 6 )=(-l .7 ,3.0) for 
cluster stars, which translates into a heliocentric tangential ve- 
locity of -300 km s _1 at 7.5 kpc. This is somewhat high, but an 
acceptable value. The observed behavior of proper motions in 
the field of VVV CL002 can therefore reflect the presence of a 
comoving group of stars. 

3.8. The nature of the cluster candidates 

Although the observed radial density profiles are well- 
-represented by the King law, which is typical of stellar clus- 
ters, and the stellar overdensity associated with VVV CL002 and 
VVV CL003 is statistically highly significant, we will not rely 
on stellar counts alone to conclude about the nature of the three 
cluster candidate s, because they can sometimes be misleading 
(see, for example, [Villanova et al.|2004||Moni Bidin et al.|2010| >. 
Unfortunately, the proper motions do not provide useful infor- 
mation about this. 

The case of VVV CL003 is straightforward, because the 
cluster sequence in the CMD is clearly distinct from the field 
one, hence this object cannot be only a fluctuation of the field 
stellar density. The redder putative cluster sequence cannot be 
ascribed to a local fluctuation of reddening, because in this case 
the extinction would be higher: a lack of stars would thus be ex- 
pected instead of the observed significantly higher stellar counts. 
The reddening-free CMD also supports this conclusion. 

The significance of the stellar overdensity detected in the 
direction of VVV CL002 is the highest among the three can- 
didates. As an additional proof that the cluster and field stars 
originate from different populations, we derived the spatial den- 
sity distribution of the stars with K s <11 located within \A(J - 
K s )\ <0.07 from the isochrone plotted in Figure[l2| which fits the 
deco ntaminated cluster sequence with the parameters derived in 
g3.9| The resulting density map is shown in the upper panel of 
Figure [T0| where each 2" x 2" pixel is assigned the average den- 
sity in the 20" x 20" square centered on it. The local overeden- 
sity associated to the cluster is clearly visible at (x,y)~(250,750). 
The procedure was then repeated, shifting the same isochrone 



by 0.14 magnitudes toward the red. The map is shown in the 
lower panel of the same figure, where the cluster overdensity is 
barely visible, lost in the patchy structures of the field. A ran- 
dom overdensity of stars should be observed in any region of 
the CMD where the field is found though. This demonstrates 
that VVV CL002 is not a simple fluctuation of the field density. 
Still, we could be in the presence of a narrow, deep window in 
the interstellar clouds, where a higher quantity of stars, on av- 
erage bluer because less reddened, would be found. However, 
the presence of a similar color offset even in the reddening-free 
CMD argues against this hypothesis. Therefore, VVV CL002 
and VVV CL003 are most likely two new stellar clusters. 

Even VVV CL004 is a simple case, because all evidence 
points at a random fluctuation of field stars. Indeed, local stel- 
lar clustering is not very significant, its stars are photometrically 
indistinguishable from the field population, and its radial density 
profile, although well fitted with a King profile, is quite peculiar. 
Nevertheless, VVV CL004 could still be an open cluster rem- 
nant (OCR), a dynamically evolved cluster in the last stage be- 
fore dissolution. Open cluster remnants should appear as small 
group s of bright stars, depleted of the fa inter main-sequence ob- 
jects ( |de La Fuente Marcos|1997[[l998| ), as can be observed in 
the field-decontaminated CMD of VVV CL004. 

In the following section we will estimate the physical pa- 
rameters of VVV CL002 and VVV CL003. The results are 
summarized in Table [3] A similar analysis is not presented for 
VVV CL004, both because the object is most probably not a real 
cluster, and because the stars in the decontaminated CMD are too 
few and dispersed for a reliable estimate. The most satisfactory 
isochrone fit of the VVV CL004 CMD is shown in Figure [12] 
It was obtained assuming [Fe/H] = -1.1, [ar/Fe] = +0.30, 
(m-M)* B = 15.4, E(J - K s )=0.82, and an age of 7.5 Gyr. It is not 
the only solution, anyway, and other similar fits can be obtained 
with a very different set of parameters. 



The transformations of Grocholski & Sara jedini (2002) were 
used to link the theoretical K values to the observed ^ s magni- 
tudes. We relied on key photometric quantities measured on the 
decontaminated CMDs, through t he calibrations of |Valenti et al.| 
(|2004b| ) and |Ferraro et af] pUfe] hereafter F06). In particular, 
F06 define two sets of equations for disk/halo-like and bulge-like 
objects, depending on the assumed trend of the ^-enhancement 
with metallicity. Indeed, the a-element abundance affects the lo- 
cation and curvature of the RGB, which are better described as 
a function of the so-called global metallicity ([M/H]) which, for 
[Fe/H] < -1, can be expressed as 



[M/H] = [Fe/H] + log(0.638/ a + 0.362), 



(4) 



whe re f a is the a-elements enhancement factor (Sal aris et al 
|1993| ). Thi s relation, however, breaks down with increasing 
metallicity ( [VandenBerg et al.|2000[|Kim et al.|2002| ). 

3.9. VVV CL002 parameters 

The decontaminated CMD of VVV CL002 shows a clump of 
stars at K s »14. We interpret this as a signature of the clus- 
ter's horizontal branch red clump. Following the prescriptions 
of |Kuchinski & Frogel| fl995] ) and |Kuchinski et al.| < [T995] >, we 
linearly fitted the upper RGB from 0.6 mag above the clump 
up to the brightest stars, obtaining a slope of -0.110. The fit 
is shown in Figure [6] This value is very similar to the slope 
of the RGB of NG C6528 ( |Ferraro et aLl|2000t and Lillerl 
(Valen tTet al.|2010| ), two metal-rich bulge globular clusters en- 
riched in Qf-elements ([ar/Fe] « +0.30). NGC 6528 is the bulge 
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Fig. 10. Upper panel. Spatial density of stars with K s < 17 lo- 
cated within |A(7 - ^ s )| < 0.07 from the isochrone of Figure 12 
for VVV CL002. The coordinates are given in pixels, and the 
density is in units of stars per arcmin 2 . Lower panel: same as 
the upper panel, but with an isochrone shifted 0.14 magnitudes 
toward the red. 



GC ri chest in iron, with [Fe/H] estimates in the range -0.1- 
+0.1 ( [Zoccali et al. | [20041 JOriglia et al.||2005[ |Carretta et aL 
2009] ), while |Origlia et aLl ( |2001j ) report [Fe/H] = -0.30 for 
Liller 1 . For VVV CL002 we used the F06 equations for bulge 
systems, because of its location very close to the Galactic cen- 
ter (see later), obtaining [M/H] = -0.16 ± 0.20 and, assuming 
[of/Fe] = +0.30, [Fe/H] = -0.4 + 0.2 from Equation Q. The 
errors were derived from the uncertainties in the equation con- 
stants estimated by Fer raro et al.| ( |2000| ), which dominate the 
error budget when compared to the uncertainties in the derived 
slope. Had we used the equations for disk/halo objects, we would 
have obtained [M/H] = -0.24 and [Fe/H] = -0.45, consistent 
with the adopted values. 

The mean magnitude and color of the stars in the red clump 
of VVV CL002 are K s = 13.80 + 0.15 and (J-K s ) = 2.17 + 0.04, 
where the rms of the measurements are taken as uncertainties on 
the average values. The absolute magnitude of the red clump 
weakl y depends on the age and metallicity of the stellar popu- 
lation. |Ferraro~eTaL ( |2006| ) propose M^ =-l-40 and a conser- 
take this eff< 



vative error to take this elfect into account, but still their cal- 
ibration is more suited for intermediate metal-poor, old clus- 
ters. VVV CL002 is at the upper limit of the metallicity range 
considered by |F06[ and its red clump is therefore expected to 
be brighter, as indicated by the th eoretical valu es tabulated by 
Salaris & Girardil (|2002 ). Indeed , |Alves| ([2000 ) report a mean 
value 0.2 magnitudes brighter than F06 for their solar-metallicity 
field sample. We fin ally adopted M^q=-1.50 + 0.15 which, as 
shown by Table 1 of Salaris & Girardi| ( |2002| ), encompasses the 



theoretical expectations for stars at [Fe/H] = -0.4 in the age 
interval from 1 to 1 1 Gyr. We thus obtain the apparent distance 
modulus (m-M)^ = 15.34 + 0.21. From Equation 3.4 of [F06| 
it follows that the apparent magnitude at the unreddened color 
(/ - K s ) =0J is ^ s = 13.52. The observed color at this magni- 
tude was obtained with a linear fit of the cluster RGB, restricted 
to the range ^ s =13-14 to limit the effects of the RGB curva- 
ture. The comparison of this value with the expectation for null 
reddening gave E(J - K s ) = 1.50+0.15. The uncertainty asso- 
ciated to the fit returned a too small error (+0.03 magnitudes), 
while the reddening of Galactic GCs is usu ally uncertain to at 
least 10% (see, for example, the discussion in |Harris|1996[ 2010 
version). The error on E(J - K s ) was therefore arbitrarily fixed 
to the more realistic value given above. The reddening agrees 
well with t he value o btained by horizontally shifting a PGPUC 
isochrone ( |Valcarce|2011| Valcarce et al. 2011, in preparation), 
with the required metallicity, ^-enhancement, and apparent dis- 
tance modulus, to find the best fit of the obs erved RGB. Through 
the transformations of |Cardelli et al.| ( [T989] ), the derived redden- 
ing translates into E(B - V)=2.88+0.29, A^ s = 1.02+0.10, im- 
plying a true distance modulus (m-M)o= 14.32+0.23 and a dis- 
tance d=7. 3+0.9 kpc. Had we used the mean intrinsic colo r of 
the red clump proposed by Lopez-Corred oira et al.| ( |2002| ), we 
would have obtained E(/ - K s ) = 1.46+0.15, in good agree- 
ment with our estimate, but in this case the final isochrone fit 



of the cluster CMD is of lesser quality. Schleg el et al.| (|1998) , s 
maps give E(B - V)=3.47 in the direction of VVV CL002, but 
this higher value should not be surprising, because it refers to 
the total Galactic absorption, while the cluster should be located 
only half way to the end of the a bsorbing gas column. In any 
case, the corrections proposed by |Bonifacio et al.| ( |2000| ) give 
E(B - V)=2.29, a value incompatible with our results. 

At a distance of 7. 3 kp c, the tidal and half-light angular 
radii of VVV CL002 ( g3.3| ) translate into the physical dimen- 
sions 3.8+0.5 and 1.6+0.3 pc, respectively. From the estimated 
distance and the Galactic coordinates (l,b), assuming a solar 
Galactocentric distance R©=8 kpc, we also derived the distance 
from the Galactic center R G c=0.6+0.8 kpc, and the height above 
the Galactic plane z=l 13+12 pc. 

The cluster age was determined by fitting the lower part 
of the decontaminated CMD with a PGPUC isochrone with 
the metallicity, ^-enhancement, distance modulus, and redden- 
ing determined above. In Figure [12| we show this fit with two 
isochrones of different age. We found that a young age (~6.5 
Gyr) was required to fit the feature at (K s , (J-K s )) « (17, 2). This 
clump of stars, closely resembling the cluster's sub-giant branch 
(SGB) and upper turn-off, was always observed when varying 
the cluster area and the comparison field in the decontamination 
procedure. However, we regard this age estimate as uncertain, 
because it completely relies on stars at the fainter end of our di- 
agram, but it can be regarded as a lower limit: had the cluster 
been younger, the SGB would have been brighter, and thus more 
clearly observed in our photometry. Deeper data unveiling the 
turnoff region are required to confirm these results. It must also 
be remembered that we used isochrones with Y=0.245, but a 
different helium abundance can affect the age estimate, because 
stars with higher helium content evolve faster. 

The cluster distance and all other distance-dependent param- 
eters such as Rgc> z, and the physical size, depend on the as- 
sumed absolute magnitude of the red clump. As an independent 
test of our results, we estimated the distance mo dulus differen- 
tially with respect to 47Tucanae ([Fe/H] = -0.7, Gratton et al. 
200 3), and NGC 6528 ([ Fe/H] « 0), assuming the photometry 
by |Valenti et al.| ([2004a) to derive the shift in color and mag- 
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nitude between them and VVV CL002, and the |Cardelli et all 
(1989) extinction coefficients. Despite the different metallicity 
of the two comparison clusters, the resulting distance modulus 
differs by only one tenth of magnitude in the two cases. We ob- 
tain (m-M) =14.16, and a distance d=6.8 kpc. that agrees within 
errors with our estimate. The difference is c aused by the old age 
of the two comparis on clusters (10-11 Gyr, |Gratton et al.|2003| 
|Ortolani et al.|200l| : the distance modulus derived with this em- 
pirical metho d is 0. 16 magnitudes lowe r than what we previously 
obtained, and Salaris & Girardi (2002 ) predict that the red clump 
of a population of similar age and [Fe/H] = -0.4 would be 0.12- 
0.15 magnitudes fainter than our assumption. We thus conclude 
that, although the two estimates agree well, a distance shorter by 
about 0.5 kpc should be preferred if VVV CL002 is an old GC. 

We estimated the total luminosity of the cluster integrating 
the flux collected within increasing apertures up to the tidal ra- 
dius, and subtracting the expected field cont ribution. The mea- 
surements were performed with SExtractor ( |Bertin & Arnouts] 
1996 ) on the stacked tile in the /-band, because many satu- 
rated stars prevented a reliable estimate in the ^ s -band frame. 
We thus found a total magnitude J tot = 12.5+0.2, and from the 
mean cluster color (/ - K2) = 2.3+0.1 and the apparent dis- 
tance modulus, we obtain M K s = -5.1+0.3, which translates 
into My = -3.4+0.3 assuming the synthetic colors of an old 



stellar population calculated by Leith erer et al.| ( |1999] ) 



3.10. VVV CL003 parameters 

The interpretation of the decontaminated CMD of VVV CL003 
must take into account the problems that its photometry pre- 
sented. The diagram reveals a clear sequence of stars distinct 
from the field population, resembling a cluster RGB. A sparse 
group of highly reddened stars is also visible to the right of 
this sequence. Many of them overlap with the cluster RGB in 
the reddening-free CMD and are probably the result of dif- 
ferential reddening, while some are most probably foreground 
cool dwarfs, because they still keep their redder color even in 
Figure [8] On the bluer side, a residual contamination is found. 
There is a lack of stars fainter than K s =16, probably because al- 
most the entire cluster area falls in the CCD border, where the 
data were collected in only one of the two jittered exposures 
forming a frame. The cluster RGB and the redder objects are 
visible for any definition of cluster and field area, provided that 
the latter is not more distant than 2-3' from the center, while the 
residual contamination may change in color and quantity of stars 
depending on the location of the comparison field. 

A dense group of stars in the cluster RGB is found at 7^=15. 
This feature is more easily observed plotting the histogram of the 
distribution of stars along the sequence, as shown in Figure [TT] 
In this plot we considered only stars within 0.2 magnitudes in 
color from the RGB ridge line, defined as the isochrone plotted 
in Figure [12] Unfortunately, the field and cluster RGBs are close 
in this part of the CMD, and the behavior of this feature (and 
even its detection) strongly depends on the choice of the com- 
parison field. However, an enhanced residual contamination ap- 
peared at the same magnitude and a bluer color when the clump 
in the RGB was not found, indicating that the field RGB was 
too red, partially overlapping and erasing the cluster clump. We 
will therefore assume this group of stars is the cluster horizontal 
branch red clump, but we caution that the subsequent analysis is 
based on this detection, which is not clear beyond any doubt. 

Repeating the same procedures used for VVV CL002 ( g3.9| ), 
the slope of the upper RGB estimated with the fit shown in the 
right panel of Figure [6] is -0.117. This value is lower than the 
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Fig. 11. Luminosity function of RGB stars in VVV CL003. 



slope of a ny of the 27 Galactic GCs studied by Ferrar o et al. 
(2000) and |Valenti et al.|p004b"l|2010| , indicating a high metal 
content. The location of VVV CL003 in the Galactic disk outside 



the bulge (see below) led us to adopt the F06 calibration for disk 
objects and no a-enhancement, obtaining [M/H] = [Fe/H] = 
-0. 1 +0.2. For the absolute magnitude of the red clump we opted 
for the value proposed by Alves (2000), given the high metallic- 
ity derived from the RGB slope. We thus obtained an apparent 
distance modulus (m-M)*. =16.58+0. 12, E(J-K S ) = 1.48+0.15 
from Equatio n 2.4 of [FO^j and E(B - V)=2.85 +0.29. Using the 
calibration of |L6pez-Corredoira et al.| ( |2002|), we would have 
obtained E(/ - K S )=1A6. On the contrary ^Schlegel et al.| ( [T998] ) 
give E(B - V)=3.34 in the direction of VVV CL003. As in the 
case of VVV CL002, this value i s too high compared w ith our 
estimates, while the corrections of Bonifaci o et al.| ( [2000| ) reduce 
it too much, returning E(B - V)=2.19. However, the resolution 
of the Schlegel et al.| ( |1998| ) maps is too low for this region of 
the sky of strong reddening gradients, and it is not surprising 
that their average value in the resolution element does not coin- 
cide with the local absorption. It could also be that the relation 
R=Ay/E(i? - V)=3.1, required to transform the reddening in dif- 
ferent bands through the relations of |Cardelli et al.| ([T989 ), is not 
valid in this highly obscured part of the sky. 

The reddening in the infrared bands implies A^ s = 
1.01+0.10, a true distance modulus (m-M) = 15.57+0. 16, and 
a distance d=13+1.0 kpc. With simple geometrical arguments 
based on the (l,b) coordinates given in Table [T] we obtain the 
Galactic height z=166+14 pc, and the Galactocentric distance 
Rgc=5. 0+1.0 kpc. At a heliocentric distance of 13 kpc, the an- 
gular size of VVV CL003 translates into the physical quantities 
r t = 6.8+0.4 pc and r h = 2.3+0.4 pc. 

An age estimate for VVV CL003 is impossible, because the 
CMD does not reach the SGB. We found that the observed RGB 
could be satisfactorily fitted with isochrones of virtually any age, 
although lower values required a progressively higher reddening 
to match the observed sequence. For example, an isochrone of 
2 Gyr better fitted the cluster CMD if E(J - K s )=1.56 is as- 
sumed. This is not necessarily a reason to prefer an older age 
for V VV CL003, because it could simply reflect the fact that the 
|F06| relation used to estimate the reddening was calibrated on old 
stellar populations. 

The high inters tella r reddening gradient and the problematic 
photometric data ( g2.2| ) prevented a reliable estimate of the clus- 
ter total luminosity. 



4. Discussion and conclusions 

The age of VVV CL002 is uncertain, but the proposed lower 
limit of 6.5 Gyr, its subsolar metallicity, and its height above the 
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Table 3. Derived parameters of the two cluster candidates. 





VVV CL002 


VVV CL003 




0.75+0.10 


0.6+0.1 


r h (pc) 


1.6+0.3 


2.3+0.4 


r t (') 


1.8+0.1 


1.8+0.1 


r t (pc) 


3.8+0.5 


6.8+0.4 


c = log (r h /r c ) 


0.65+0.26 


0.56+0.21 


(m-M)o 


14.32+0.23 


15.57+0.16 


d (kpc) 


7.3+0.9 


13.0+1.0 


R GC (kpc) 


0.7+0.9 


5.0+1.0 


z(pc) 


113+13 


166+14 


E(7 - K s ) 


1.50+0.15 


1.48+0.15 


E(B - V) 


2.88+0.29 


2.85 + 0.29 


[M/H] 


-0.16 + 0.2 


-0.1+0.2 


[Fe/H] 


-0.4 + 0.2 


-0.1 + 0.2 


My 


-3.4+0.3 




Age (Gyr) 


>6.5 




VVV CL002 


VVV CL003 


VVV 




Fig. 12. Decontaminated CMDs of the three clusters, with the 
best-fitting PGPUC isochrones defined in the text overplot- 
ted. The zero-age horizontal branch (ZAHB) in all panels of 
the figure was calculated with the helium abundance Y=0.245. 
For VVV CL002, two isochrones at 6.5 Gyr (thick curve) and 
8.5 Gyr (dashed curve) are shown. 



Gala ctic plane, twice the scale height of young objects ( |Joshi| 
2005 ), argue in favor of an old cluster. On the other hand, |FrieT| 
( [1995] ) showed that old open clusters are completely absent in 
the inner ~7 kpc from the Galactic center. Because of its location 
very close to the Galactic center, VVV CL002 is therefore most 
likely a GC. 

So far, HP 1 is the closest kno wn GC to the Galactic center 
=0.5 kpc |Barbuy et al.|2006 



(Rgc 



and only six GCs have been 



discovered in the inner 1 kpc (Harris 1996). VVV CL002, with 
an estimated Galactocentric distance of 0.7 kpc, is therefore one 
of the innermost Galactic GCs. Curiously, the density profile of 
VVV CL002 suggests the presence of tidal tails, a behavior more 



easily observed whe n the cluster is close to the perigalacticon 
(Kiipp eret al.|2010a) . If the connection between tidal tails and 
perigalacticon were to be confirmed even for objects subject to 
the intense stress of the inner Galaxy, then the cluster orbit would 
be confined to the very central region of the Galaxy. 

VVV CL002 has a global metallicity close to NGC 6528 
and Liller 1, as indicated by the slope of the upper RGB. These 
two objects are also bulge GCs ver y close to the Galactic center 
(Rgc=0-6 and 0.8 kpc, respectively; Valent fet al . |20 1 1 [Feltzing 



|& Johnson||2002t and, as in the case of VVV CL002, their 
stars also are barely dis tinguished fro m the background bulge 
population in the CMD (Davidge 2000). However, the similar- 
ity cannot be drawn further, because their appearance is very 
distinct. In Figure 13 we show a visual comparison between 
the three objects, with 7<r s -band images from the VSA archive: 
VVV CL002 is clearly sparser and less massive, resembling a 
Palomar-like object more than a highly concentrated bulge GC 
such as NGC 6528 and Liller 1 . Indeed, our estimate of its to- 
tal luminosity is -3.5 magnitudes fainter than these two objects 
(|Webbink||1985[|Peterson & Reed||T987l |van den Bergh et al 



[T99T}|Mallen-Ornelas & Djorgov ski[T993l ). VVV CL002 is more 



affine to moderately metal-rich, low-mass, and sparse GCs such 
as Whiting 1 and Palomar 1, whose ph ysical size are compara 
ble (r h =1.9 and 1.5-2.2 pc, respectively ; |Harris| 1 996| [Ro senberg 



etal, 



et al. 



1998 ). These two objects are relatively young (Rosenberg 



1998; Carraro 2005]) , and are only slight ly more metal- 



poor ( Carraro et al.|2007b||Monaco et al.||2011| ), although they 



are one magnitude brighter (Harris 



1996|). Palomar 12 also has 



19891 |Geisler"eTaLl[2007l ) 



a similar young age ( [Stetson et al. 

but has a lower me tallicity ([Cohen 2004]), and it is physically 
larger (Harris 1996). It is worth noting that both Whiting 1 and 
Palomar 12 are associated with the Sagittarius dwarf spheroidal 
galax y ( |Irwin|[T999} |Carraro et al.||2007bl |Law & Majewskij 
|2010| ), while the extragalactic origin of Palomar 1 has been re- 
peate dly proposed although never proven (e.g., |Crane et al.| 
|200^|Frinchaboy et al.|2004[|Belokurov et al.|2007| )~ 

VVV CL003 is surely a physical cluster of stars, located 
in the Galactic disk beyond the center, outside the bulge, at a 
Galactocentric distance of 5 kpc. Its high metallicity and its posi- 
tion in the Galaxy are more typical of an open cluster, although a 
GC cannot be excluded. Deeper photometry, reaching at least the 
cluster SGB two magnitudes fainter than the limit reached with 
our data, is needed to estimate its age and, thus, better clarify 
its nature. Interestingly enough, to our knowledge it is the first 
stellar cluster ever discovered on the "dark side of the Galaxy", 
i.e. the Galactic disk on the other side of the c enter. The Open 
Cluster Databasf](WEBDA, (Mermilhod|l 996 ) has no entry in 
the direction of the bulge more distant than a few kpc. 

VVV CL004 is most probably only an asterism, and as such 
it deserves no further attention. Nevertheless, the possibility of 
an OCR remains open, and this kind of object is of great inter- 
est, because we are still lacking observational evidence of the 
last stages of the cluster dynamical evolution. Indeed, OCRs are 
very elusive and easily confused with rand om field fluctuations 
and, despite a lar ge list of candidates (e.g., |Pavani & Bica|2007[ 
Bica et al.|2001| ), only one genuine OCR has been to date spec- 



troscopically confirmed (NGC 1901, |Carraro et al.|20 07a). 
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Fig. 13. 2f4 x 2f4 logarithmic gray-scale ^ s -band images of (from left to right) NGC 6528, Liller 1 , and VVV CL002, from the VVV 
archive. North is up and East is to the right. 
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